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We uncovered a novel role for the spliceosome in regulating mRNA expression levels that involves splicing 
coupled to RNA decay, which we refer to as spliceosome-mediated decay (SMD). Our transcriptome-wide studies 
identified numerous transcripts that are not known to have introns but are spliced by the spliceosome at canonical 
splice sites in Saccharomyces cerevisiae. Products of SMD are primarily degraded by the nuclear RNA 
surveillance machinery. We demonstrate that SMD can significantly down-regulate mRNA levels; splicing at 
canonical splice sites in the bromodomain factor 2 [BDF2] transcript reduced transcript levels roughly threefold by 
generating unstable products that are rapidly degraded by the nuclear surveillance machinery. Regulation of BDF2 
mRNA levels by SMD requires Bdf 1, a functionally redundant Bdf2 paralog that plays a role in recruiting the 
spliceosome to the BDF2 mRNA. Interestingly, mutating BDF2 5' splice site and branch point consensus 
sequences partially suppresses the bdflA temperature-sensitive phenotype, suggesting that maintaining proper 
levels of Bdf2 via SMD is biologically important. We propose that the spliceosome can also repress protein-coding 
gene expression by promoting nuclear turnover of spliced RNA products and provide an insight for coordinated 
regulation of Bdf 1 and Bdf2 levels in the cell. 
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The biogenesis of pre-mRNA involves multiple process- 
ing reactions, including 5 ' capping, splicing, cleavage, and 
polyadenylation at the 3' end, leading to production of 
translationally competent mRNA (Moore and Proudfoot 
2009). In nuclear pre-mRNA splicing, the excision of 
introns is catalyzed by the spliceosome, a ribonucleo- 
protein machine comprising five small nuclear RNAs 
(snRNAs) and >100 conserved proteins (Wahl et al. 2009). 
Spliceosomal snRNAs contain a conserved PuAU 4 _ 6 GPu 
sequence called the Sm site, which provides a platform 
for the assembly of the heteroheptameric Sm complex, 
comprised of seven Sm proteins: SmB/B', SmDl, SmD2, 
SmD3, SmE, SmF, and SmG (Beggs 2005). These form 
a ring, referred to as the Sm complex, around the Sm site 
of the spliceosomal Ul, U2, U4, and U5 snRNAs, and play 
roles in multiple aspects of small nuclear ribonucleopro- 
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tein (snRNP) biogenesis. The U6 snRNA associates with 
a structurally related set of seven Lsm (like Sm) proteins. 
Being a part of several snRNPs, the Sm complex was 
shown to facilitate assembly of the spliceosome on pre- 
mRNA (Zhang et al. 2001) and play a role in multiple 
aspects of snRNA biogenesis, such as cellular localiza- 
tion, processing, and stability. 

The spliceosomal snRNPs and multiple non-snRNP 
proteins assemble ^transcriptionally on pre-mRNAs 
through recognition of the 5' splice site (5'ss), the branch 
point (BP), and the 3'ss to form the spliceosome (Carrillo 
Oesterreich et al. 2011). This involves interactions of the 
spliceosomal components with the 5' cap-binding com- 
plex and the C-terminal domain (CTD) of RNA poly- 
merase II (RNA Pol II). Splicing factors such as yeast 
Prp40 and several human SR proteins were implicated in 
mediating interactions with the phosphorylated CTD 
(Morris and Greenleaf 2000; Kotovic et al. 2003). Pre- 
vious analyses of spliceosomal recruitment to RNA Pol 
II-transcribed genes were performed by studying the 
genome-wide distribution of the individual components 
of the spliceosome using chromatin immunoprecipitation 
(ChIP) (Kotovic et al. 2003; Moore et al. 2006; Tardiff et al. 
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2006). Surprisingly, these studies revealed that spliceosomal 
components are also recruited to some genes whose 
transcripts are not known to be spliced. It was recently 
shown that several nonintronic transcripts can be spliced; 
however, it remained unclear whether this had any 
regulatory role or was the result of stochastic splicing 
events that occur due to the lack of specificity in the 
spliceosome recruitment (Harigaya and Parker 2012). 

Unexpectedly, by sequencing SmDl -associated RNAs, 
we identified many RNA Pol II transcripts that are not 
known to be spliced, indicating that spliceosome assem- 
bly on nonintronic mRNAs may be a common phenom- 
enon. To understand the functional significance of this 
discovery, we analyzed the effect of splicing mutants and 
strains defective in RNA surveillance on the protein- 
coding transcriptome. Remarkably, these mutants dis- 
played a marked up -regulation of a number of intronless 
mRNAs, including bromodomain factor 2 (BDF2) mRNA, 
which encodes a bromodomain factor, and OYE3 mRNA, 
which encodes an NADPH oxidoreductase. Interestingly, 
it was previously shown that BDF2 overexpression is toxic, 
implying that it could be important to down-regulate levels 
of this transcript in the cell (Yoshikawa et al. 2011; Fu 
et al. 2013). We demonstrate that splicing of these mRNAs 
at consensus splice sites generates unstable products that 
are primarily degraded by the nuclear RNA surveillance 
machinery. This spliceo some-mediated decay (hereafter 
referred to as SMD) of BDF2 mRNA is dependent on Bdf 1, 
another bromodomain-containing protein, providing a 
plausible explanation of how expression of these paral- 
ogous genes is regulated. Mutating BDF2 5'ss and BP 
consensus sequences suppressed bdflA growth pheno- 
types, suggesting that maintaining proper levels of Bdf2 
via SMD is biologically important. We propose a model in 
which SMD regulates coordinated expression of Bdf 1 and 
Bdf 2. Collectively, our results have revealed a new role for 
the spliceosome in the regulation of mRNA expression. 

Results 

Purification and analyses of the Sm complex-bound 
RNAs 

To identify RNA targets of the spliceosome, we purified 
spliceosome-associated RNAs by affinity chromatogra- 
phy using a strain expressing tandem affinity purification 
(TAP)-tagged SmDl, a protein required for stabilizing the 
Ul snRNA-pre-mRNA interaction (Zhang et al. 2001 ). As 
a negative control, we used a nontagged parental strain. 
SmDl-coprecipitated RNA was extracted and analyzed 
by Northern blotting and high-throughput RNA sequenc- 
ing (RNA-seq). Methylene blue staining of total RNA 
showed a significant depletion of ribosomal RNAs (rRNAs) 
in the SmDl pull-down compared with input (Supplemen- 
tal Fig. 1A). Conversely, Northern analysis revealed sub- 
stantial enrichment of the Sm site-containing RNAs Ul 
(Supplemental Fig. 1A) and telomerase (TLC1) (data not 
shown), validating the biological relevance of this ap- 
proach. Ul was undetectable in the mock samples, dem- 
onstrating the specificity of the results (Supplemental Fig. 



1A, lanes 5,7). Mass spectrometric analysis of coprecipitated 
proteins identified other Sm complex members (SmB, 
SmDl, SmD2, SmE, and SmG) and many spliceosome 
components (Supplemental Fig. IB). 

Bioinformatics analysis of the RNA-seq data (Supplemen- 
tal Table 1) revealed that 35%-40% of uniquely mapped 
cDNAs corresponded to snRNAs directly bound by the 
Sm complex, with the highly abundant U2 snRNA being 
the most highly represented RNA. U6 snRNA was also 
present at low levels, indicating that the U4/U6 particle 
as well as late spliceosome assembly intermediates 
also coprecipitated (Will and Luhrmann 2011). Consis- 
tent with our Northern blotting data, the low-abundant 
TLC1 RNA that contains an Sm-binding site (Seto et al. 
1999) was substantially enriched relative to other non- 
coding RNAs (ncRNAs), such as SCR1, which coprecipi- 
tated nonspecifically at low levels compared with control 
samples (Supplemental Table 1 ). Notably, the snR190 box 
C/D small nucleolar RNA (snoRNA) was also enriched, 
consistent with previously published data (Zagorski et al. 
1988). Collectively, these data demonstrate the specific- 
ity and sensitivity of our analysis. 

Of the top 500 enriched mRNAs, 412 were reproduc- 
ibly identified in both data sets (Fig. 1 A). These included 
198 intron-containing mRNAs, representing —60% of the 
known spliced mRNAs in Sacchawmyces cerevisiae 
(341) (Fig. IB). It is possible that the fraction correspond- 
ing to the remaining 40% included pre-mRNAs that are 
rapidly spliced and therefore underrepresented in the 
sequencing data. Strikingly, more than half of the 412 
mRNAs identified were annotated as nonintronic (Fig. IB). 
We envisage two possible explanations for this observa- 
tion. First, the Sm complex could assemble directly on 
mRNAs independently of the spliceosome. Indeed, in 
silico analyses revealed 2558 potential Sm-binding sites 
in intronless genes and 94 in the 213 Sm-associated 
mRNAs (data not shown). Alternatively, this interaction 
could be mediated via the spliceosome, which is known 
to be recruited to a subset of intronless mRNAs (Koto vie 
et al. 2003; Tardiff et al. 2006). The relevance of these inter- 
actions, however, was not addressed, begging the question 
of whether these mRNAs are spliced by the spliceosome. 
Our in silico analyses revealed that 44 of the 213 intron- 
less mRNAs contained canonical 5'ss (GUA[U/C/A]GU) 
and BP (ACUAAC[G/A/U]| sequences in the correct 
orientation (Fig. 1C). Three were previously described 
[FBA1, UTH1, and PGI1 ) (Kotovic et al. 2003), and six had 
a putative Sm site ([A/G]AU 4/6 G[A/G]) within 50 base 
pairs (bp) upstream of the 5'ss [KAR2, KEM1, PSK1, STE6, 
and YPR045C). To validate our findings, we initially 
focused on the BDF2 mRNA, as it was most likely to be 
targeted by the spliceosome for the following reasons: (1) 
BDF2 RNA was highly enriched in both Sm-IP data sets 
comparable with intron-containing mRNAs (Supplemen- 
tal Table 1) and has predicted 5'ss, BP, and potential 3'ss 
(Zhang et al. 2007) in the correct orientation (Fig. IE). (2) 
Ul, U2, and U5 components were found to chromatin- 
immunoprecipitate over the BDF2 gene at levels compa- 
rable with a typical intron-containing gene, RPL28A (Fig. 
ID; Tardiff et al. 2006). (3) Intron tiling array data from the 
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Figure 1. Sm proteins associate with nonintronic mRNAs. [A] Venn diagram showing the overlap between the top 500 mRNAs 
enriched in two separate SmDl-TAP experiments. [B] Venn diagram showing that of the 412 genes that reproducibly coprecipitated 
with the Sm complex, 198 contained introns, and 213 were annotated as intronless genes. The pink circle shows the 412 common genes 
from Sm-TAP immunoprecipitation; of these, 213 genes are nonintronic. The green circle shows 341 annotated intronic genes in S. 
cerevisiae-, of these, 198 were found in Sm-TAP immunoprecipitation. (C) Venn diagram showing that out of all intronless genes, 835 
are predicted to contain BP and 5'ss sequences in the correct order and have a potential to recruit the spliceosome. Forty-four genes 
with the splice site consensus were present in Sm-IP. [D,E) University of California at Santa Cruz Genome Browser outputs displaying 
the Sm RNA-seq results for RPL28A, BDF2, and neighboring genes. Included are the positions of the predicted 5'ss and BP splice 
signals. Spliceosome ChIP tiling array data (Ul, U2, and U5) (Tardiff et al. 2006) and dbrIA tiling array data are indicated as "Arrays/ 7 
and log 2 changes between mutant and wild type are shown (Tardiff et al. 2006; ENCODE Project Consortium et al. 2007). 
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dbrlA strain, where the lariat intermediate is stabilized 
due to the loss of debranching enzyme, have indicated the 
presence of a putative intron in BDF2 mRNA (Zhang 
et al. 2007). 

The spliceosome endonucleolytically cleaves BDF2 
mRNA 

To determine whether expression of the intronless BDF2 
mRNA is regulated by the spliceosome, we analyzed its 
steady-state levels in the prp40-l temperature-sensitive 
strain. Prp40 is a subunit of the Ul snRNP involved in 
cotranscriptional recruitment of early splicing factors to 
pre-mRNAs and later steps of spliceosome assembly. 
Previous analysis has demonstrated global down-regula- 
tion of mRNA levels in this mutant as a result of splicing 
defects (Albulescu et al. 2012). As expected, at the non- 
permissive temperature, mRNA levels of the intron- 
containing gene RPL28A were substantially reduced 
(Fig. 2A). In contrast, in this strain, BDF2 mRNA levels 
increased roughly threefold (Fig. 2B). Similarly, BDF2 
mRNA levels were also elevated in Ul -depleted cells 
and prp2-l (Fig. 2C, lanes 1-4) andprp42-l (data not shown) 
mutants, where the first catalytic step of the splicing 
reaction is affected. A mutant defective in catalysis of the 
second transesterification reaction (prpl7-l) (Jones et al. 
1995; Noble and Guthrie 1996; Pleiss et al. 2007) also 
accumulated BDF2 mRNA (Fig. 2C, lanes 5,6). Finally, 
a 2-base substitution introduced into the predicted 5'ss 
(Fig. 2D) also increased BDF2 mRNA levels (Fig. 2E). We 
conclude that the spliceosome plays a direct role in down- 
regulating BDF2 expression. To substantiate our results, 
we also performed Northern blots to measure levels 
of splicing intermediates in strains defective in splic- 
ing (prp40-l), nuclear RNA surveillance (rrp6A), and/or 



debranching of the intron lariat (dbrlA). Using a probe 
complementary to a region downstream from the pre- 
dicted 5' exon (probe 2) (Fig. 3 A), we could readily detect a 
lariat in the dbrlA strain (Fig. 3B [lane 5], C [lane 2]). Cells 
lacking both Dbrl and Rrp6 also accumulated a 3' ex- 
tended lariat intermediate (Fig. 3C [lanes 3,6], D [lane 3]). 
Excised exon 1, a product of the first splicing reaction, 
could also be detected in the rrp6A strain (Fig. 3E, lane 4), 
an exonuclease-defective dis3 mutant, and a ratl-1 tem- 
perature-sensitive mutant (Fig. 3E, lanes 6,8). We note 
that exon 1 fragments were heterogeneous in length, 
which could be the result of exosome-mediated degrada- 
tion and/or low-fidelity splicing events. Deletion of the 
cytoplasmic exonuclease Xrnl and the nonsense-medi- 
ated decay (NMD) factor Upf 1 had no detectable effect on 
the stability of the exon 1 splicing products (Fig. 3E, lanes 
5,7). Altogether, our data indicate that spliceosomal BDF2 
cleavage products are mainly degraded in the nucleus by 
both 5'-3' (Rati), and 3'-5' (Rrp6 and Dis3) exonucleases. 
In contrast to BDF2, we could not detect accumulation of 
spliced RPL28A 5'exons (Fig. 2 A), which is consistent with 
published observations for typical intronic genes (Egecioglu 
and Chanfreau 2011). This is most likely due to the fast 
kinetics of the second splicing reaction. 

Accumulation of BDF2 5' exon was also dependent on 
intact 5'ss and the BP sequences (Fig. 3F, lanes 3-6). In 
contrast, mutations in three predicted 3'ss had no notice- 
able effect on exon 1 levels or BDF2 mRNA levels in 
general (Supplemental Fig. 2 A; data not shown), indicat- 
ing that splicing at these sites does not occur or is very 
inefficient. We conclude that splicing products generated 
using these splice elements are products of a single 
splicing step and are degraded by the nuclear surveillance 
system. This indicates that the second step of splicing in 
the BDF2 mRNA is inefficient, as it was reported for the 
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Figure 2. BDF2 mRNA levels are up-regulated in 
the splicing mutants or upon 5'ss mutation. [A,B) 
Northern blot analysis on total RNA isolated from 
the parental strain (wild type [wt]) and pip40-l, 
prp40-l rrp6A, and rrp6A strains. RNA was ex- 
tracted 2 h after the shift to 37°C, resolved on 
agarose [A) or polyacrylamide gels [B), transferred 
to positively charged membrane, and probed for 
RPL28A (using 5' exon-specific probe) and BDF2 
mRNAs, respectively. The positions of pre-mRNA, 
mature mRNA, and 5' exon (for RPL28A) species 
are indicated. ADH1 (A) or methylene blue-stained 
25S and 18S rRNAs (B) are shown as loading con- 
trols. (C) Splicing mutants accumulate BDF2 mRNA. 
To analyze BDF2 mRNA levels, total RNA was 
extracted from Ul snRNA-depleted cells (lane 2) 
and prp2-l (lanes 3,4) and prpl7-l (lanes 5,6) mu- 
tants. Cultures were shifted to nonpermissive tem- 
perature (37°C) for 1 h (prp2-l and pip 17-1) or grown 
in the presence of glucose for 8 h to allow for Ul 
depletion. (D) Schematic representation of the pre- 
dicted splicing signals in BDF2. The purple nucleo- 
tides indicate the substitutions that were made in 
the 5'ss (red) to generate the 5'ss* construct. (£) 
Northern blot analysis of BDF2 mRNA levels in 
wild-type [BDF2) or 5'ss mutant (5'ss*) cells. 
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Figure 3. The spliceosome performs a one-step 
splicing reaction on BDF2 mRNA. (A) Schematic 
representation of the BDF2 gene. The approximate 
positions of probes used to detect BDF2 mRNA and 
spliced products are indicated. (B,C) A BDF2 intron- 
lariat accumulates in the exosome and debranching 
mutants. Northern blot analysis of the BDF2 tran- 
script was performed on RNA extracted from wild- 
type (wt), prp40-l, prp40-l np6A, np6A, dbrlA, and 
dbrlArrp6A strains. Probes 2 and 3 were used to 
detect the BDF2 mRNA and lariat-intron splicing 
intermediates. rRNAs are shown as a loading con- 
trol. A schematic representation of splicing inter- 
mediates is indicated at the side of each panel. (D) 
Primer extension analysis of the BDF2 splicing in- 
termediate in YF336, dbrlA, rrp6A dbrlA, and rrp6A 
strains. The position of the primer used in this 
experiment is indicated in the diagram. Primer 
extension products were separated on a 10% urea- 
PAGE. The arrow indicates the position of the 
primer extension product corresponding to the 3' 
extended intron-lariat. Primer extension was per- 
formed on PMA1 to control for the RNA levels. [E) 
Released BDF2 exon 1 is degraded in the nucleus. 
Northern blotting was performed on RNA extracted 
from the wild-type, prp40-l, prp40-l rrp6A, rrp6A, 
upflA, xrnlA, dis3exo-, and ratl-1 strains as de- 
scribed above. Probe 1 was used to detect the exon 1 
fragment. [F) Mutations in the 5'ss and BP abolish 
the first step of BDF2 splicing. Shown is a schematic 
representation of the mutations introduced to BDF2 
splice sites and Northern blot analysis of BDF2 
exon 1 in the splice site mutants. (G) The first step 
of splicing is affected in the prpl7-l mutant. Exon 1 
levels were analyzed by Northern blot. RNA was 
extracted from the wild-type (DIS3), prpl7-l, and 
pipll-1 dis3 exo- strains. 



telomerase RNA (Box et al. 2008; Kannan et al. 2013). BDF2 
mRNA levels also increased in the prpl7-l mutant, sug- 
gesting that BDF2 can be fully spliced (Fig. 2C); however, we 
also noticed a reduction in exon 1 levels in the pipll-1 
strain. Consistent with previous reports, this suggests that 
first-step splicing is also delayed in this mutant (Sapra 
et al. 2008), making it difficult to determine what fraction 
of BDF2 is fully spliced. 

Shortening BP-3'ss distance partially restores 
the second step of splicing in BDF2 

The distance between the BP and 3'ss was reported to be 
a factor affecting efficient 3'ss selection during the second 
catalytic step of splicing (Luukkonen and Seraphin 1997; 
Box et al. 2008; Meyer et al. 2011; Kannan et al. 2013). 
Hence, we reasoned that the unusual long distance 
between these two elements in the BDF2 mRNA might 
be responsible for the accumulation of the RNA species 
that are produced after the first step of splicing. To test 
whether BDF2 mRNA has a potential to be fully spliced, 



we generated several BDF2 constructs in which the 
distance between the BP and previously predicted 3'ss 
sites (Zhang et al. 2007) was shortened (Supplemental Fig. 
2B, A1-A3). Most of the tested constructs failed to allow 
complete splicing at detectable levels even using a sensi- 
tive technique such as RT-PCR (Supplemental Fig. 2C, 
lanes 1,3,4,6). Intriguingly, low levels of spliced BDF2 A2 
product could be detected in an Rrp6 deletion strain 
(Supplemental Fig. 2C, lane 5), indicating that this product 
is degraded by the nuclear exosome. Sequencing of the 
PCR product revealed that this indeed represented a spliced 
BDF2 mRNA in which the distal 3'ss was used (Supple- 
mental Fig. 2D). We conclude that the previously proposed 
splicing signals are suboptimally positioned in BDF2 
mRNA to promote an efficient second step of splicing. 

Both a single splicing event and complete splicing 
contribute to regulation of BDF2 levels via SMD 

While this manuscript was in preparation, the Parker 
laboratory (Harigaya and Parker 2012) reported splicing of 
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BDF2 in cells lacking Xrnl and a cytoplasmic decapping 
factor, Dcp2. However, the biological relevance of this 
splicing event had not been addressed (Harigaya and 
Parker 2012). Sequencing of the spliced product revealed 
that the 5'ss reported here and a more downstream 3'ss at 
nucleotide +1672 were used (Fig. 4A ; Harigaya and Parker 
2012). A BP consensus located 12 nucleotides (nt) upstream 
of this 3'ss was proposed; however, usage of this element 
was not experimentally validated. 

We were able to detect splicing at these sites in the 
wild-type strain,- however, no significant additional accu- 
mulation of the spliced product was in the xrnl A dcp2A 
double mutant (Fig. 4B, lanes 1,4). Deleting UPF1 also did 
not noticeably stabilize the spliced product (Fig. 4B, lane 
3), indicating that the cytoplasmic RNA surveillance 
machinery plays a minor role in the degradation of the 
spliced product. In sharp contrast, deleting RRP6 led to an 
~4. 5-fold increase in the level of the spliced product 
(Fig. 4B, lane 2). We conclude that SMD triggers the 
degradation of partially and completely spliced BDF2 
products by the nuclear surveillance machinery. 



To address the importance of this two-step splicing 
event in the regulation of BDF2 expression in more detail, 
we introduced mutations in the splice sites proposed by 
Harigaya and Parker (2012). BDF2 mRNAs with muta- 
tions in the 5'ss or the BP at +1660 (Fig. 4A) were not 
detectably spliced, indicating that these sites are indeed 
used by the spliceosome. Interestingly, although the BP 
sequence at +254 (Fig. 4A) is clearly important for the 
single-step splicing event that leads to the accumulation 
of the exon 1 fragment in the rrp6A strain (Fig. 3F), it is 
not used for the two-step splicing reaction (Fig. 4C, lane 3), 
presumably because the distance between this BP and the 
downstream 3'ss is far from optimal to support splicing 
(>1400 nt). Surprisingly, mutating the 3'ss at +1672 did not 
block splicing of BDF2, as an alternative 3'ss was used at 
+1702 (Fig. 4C, lane 5). Mutations in the 5'ss or the two 
BPs led to a twofold to threefold increase in the BDF2 
mRNA levels (Fig. 4C, lanes 2-4). Mutation in the 3'ss 
at +1672 had only a modest effect on BDF2 steady-state 
levels (Fig. 4C, lane 5), presumably because alternative 
3'ss could still be used. In contrast, the BP at +1660 and 
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Figure 4. Two-step splicing also contributes to regulation of BDF2 levels. (A) Shown is a schematic representation of the BDF2 splice 
sites, with the position relative to ATG (+1) indicated. RT product is indicated in red, and primers used for PCR are shown as black 
arrows. (5) RT-PCR results performed on oligo[d(T)]-enriched RNA extracted from the wild-type (wt), rrp6A, upflA, and xrnlA dcp2A 
strains using BDF2 and ADH1 primers. RT-PCR reactions without enzyme (RT— ) are included as a negative control. The spliced 
product is indicated. (C) RT-PCR results performed on poly(A)-enriched RNA extracted from strains expressing the wild-type [BDF2) 
and BDF2 5'ss, BP (+254), BP (+1660) (changing AAC to ACC), and 3'ss (+1672) (changing AG to AA) mutants in bdf2A up6A. (D) 
Northern blot analysis of exon 1 levels in the strains indicated. 
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the 3'ss at +1672 are not required for the single-step 
splicing event, as we could still detect accumulation of 
exon 1 in the rrp6A strain in cells expressing these BDF2 
mutations (Fig. 4D). 

We conclude that SMD of BDF2 involves (at least) two 
separate splicing events: a single cleavage event, which 
requires the BP at +254, and a two-step splicing event 
involving the BP at +1660. Both pathways use the same 5'ss, 
and spliced products generated by these pathways are 
(primarily) degraded by the nuclear surveillance machinery. 

Spliceosome recruitment to BDF2 mRNA 
is compromised in bdflA 

Bdf2 and the closely related Bdfl both bind to acetylated 
histones, albeit with different specificities (Matangkasombut 
and Buratowski 2003). These genes are genetically re- 
dundant, as only one of the two genes is necessary and 
sufficient for cell viability. Deletion of both genes is 
lethal (Matangkasombut et al. 2000). Interestingly, de- 



letion of BDF1 resulted in an approximately threefold 
increase in the levels of BDF2 mRNA (Fig. 5A ; Fu et al. 
2013), similar to what we observed in the splicing mutants 
(Fig. 2), and a twofold increase in Bdf2 protein levels (Fig. 
5B). Although Bdfl binds to a BDF2 promoter element 
(Durant and Pugh 2007; Fu et al. 2013), nuclear run-on 
experiments demonstrated that deletion of BDF1 did 
not noticeably affect BDF2 transcription (Fig. 5C), ruling 
out that Bdfl regulates Bdf2 at the transcriptional level. 
Deletion of BDF1 is known to reduce splicing of a subset 
of intron-containing transcripts (Albulescu et al. 2012), 
but the effect of BDF1 deletion on intronless mRNA 
levels was not reported. This prompted us to investigate 
whether Bdfl is required for spliceosome recruitment to 
BDF2. ChIP assays revealed a substantial Ul snRNP 
enrichment over the BDF2 gene, in agreement with 
previously published data (Fig. 5D,E [lanes 1-4], F ; Tardiff 
et al. 2006). Strikingly, deleting BDF1 dramatically re- 
duced the Ul binding to BDF2 (Fig. 5E [lanes 5-8], F). 
Similarly, deletion of BDF1 was reported to result in 
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Figure 5. Bdfl is required for the recruit- 
ment of the spliceosome at the BDF2 locus. 
(A) Northern blot analysis of BDF2 mRNA 
levels in the wild type (wt) and bdfl A. [B) 
Western blot analysis of Bdf2-TAP and Rpbl 
(3E1 antibody from Millipore) protein levels 
in the wild type and bdfl A. (C) Deletion of 
Bdfl does not noticeably affect BDF2 tran- 
scription. Transcriptional run-on experiments 
were performed to measure transcription on 
BDF2 in wild-type and bdfl A strains. RNA 
probes corresponding to BDF2 [probe 1 (exon 
1) and probe 2], ADH1 and a nontranscribed 
region on chromosome V (nc) were used for 
detection. (D) Schematic representations of 
the BDF2 gene and BDF2 mRNA with the 
predicted 5' exon (exon 1) and BR Black 
bars show the locations of the PCR products 
generated in the ChIP experiments. [E,F) 
Deletion of BDF1 reduces spliceosome re- 
cruitment to BDF2. ChIP experiments were 
performed in wild-type and bdfl A strains 
expressing TAP-tagged Yhcl (U1C), a Ul 
snRNP component. Ul enrichment was mea- 
sured by radioactive PCR using primers that 
generate the PCR fragments shown in D. PCR 
products were resolved on polyacrylamide 
gels and detected by autoradiography. Quanti- 
fication results from four independent experi- 
ments are shown in F-, the error bars indicate 
the standard error. Stars indicate position of the 
PCR product from a nontranscribed region on 
the chromosome. (G-7) Mutating BDF2 5'ss 
and BPs suppresses bdflA salt- and tempera- 
ture-sensitive phenotypes. The bdfl A strains 
expressing either BDF2 or bdf2 5'ss, 3'ss, and 
BP mutants were grown in liquid medium until 
OD(A 60 o) reached 0.5. To assay for growth, 5-jjlL 
aliquots of 10-fold serial dilutions were spotted 
onto YPD and —LEU medium in the presence 
or absence of 0.6 M NaCl. 
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compromised Ul recruitment to several intron-contain- 
ing genes, implying that Bdf 1 might play a general role in 
splicing (Albulescu et al. 2012). We conclude that BDF1 
down-regulates BDF2 mRNA levels by stimulating splic- 
ing of BDF2. 

Disrupting BDF2 splicing suppresses bdflA 
temperature- and salt- sensitive phenotypes 

To address the biological significance of regulation of 
BDF2 mRNA levels via SMD, we used a genetic approach 
that relied on growth defects observed in the bdflA strain. 
Overexpression of BDF2 was recently shown to suppress 
bdflA temperature- and salt-sensitive phenotypes in a 
dosage-dependent manner (Matangkasombut et al. 2000; 
Fu et al. 2013). We asked whether completely blocking 
splicing of BDF2 by mutating the 5'ss or blocking two- 
step splicing by mutating the BP at +1660 (see Fig. 4A) 
could at least partially suppress bdflA growth defects, 
fudging from the number of colonies in each dilution and 
the colony size, it can be concluded that at the non- 
permissive temperature (37°C), bdflA strains expressing 
BDF2 mRNAs with mutated 5'ss grew better compared 
with the bdflA cells expressing wild- type BDF2 (approx- 
imately eightfold to 10-fold) (Fig. 5G). Moreover, com- 
pletely blocking BDF2 splicing substantially enhanced 
the cells 7 resistance (~30-fold to 50-fold compared with 
bdflA BDF2) to a combination of high temperature and 
high salt concentrations (0.6 M NaCl) (Fig. 5H). Finally, 
mutating the BPs also improved growth in high salt 
conditions (Fig. 51). Collectively, these results demon- 
strate that even a relatively small change in BDF2 mRNA 
and protein levels can dramatically affect a cell's re- 
sistance to certain stress conditions, underscoring the 
impact that SMD can have on cell fitness. 

SMD may regulate expression of —1% of the intronless 
genes 

Several recent studies have described interactions be- 
tween the spliceosome and mRNAs that were not known 
to contain introns (Kotovic et al. 2003; Moore et al. 2006; 
Tardiff et al. 2006; Harigaya and Parker 2012). Our 
bioinformatics analysis indicates that >800 genes encod- 
ing intronless mRNA have canonical 5'ss and BP se- 
quences in the correct orientation (Fig. 1C), suggesting 
that the spliceosome could potentially target and regulate 
expression of hundreds of intronless mRNAs. To test this 
possibility, we performed RNA-seq analyses on rRNA- 
depleted RNA extracted from a wild-type strain and the 
temperature-sensitive prp40-l and rrp6A mutants. 

Based on our analysis of BDF2 mRNA, we reasoned 
that splicing products derived from intronless genes would 
be stabilized upon deleting Rrp6, whereas unspliced or 
full-length transcripts would accumulate in the prp40-l 
mutant. RNA was extracted from these strains 2 h after 
the shift to the nonpermissive temperature. This was 
sufficient to detect accumulation of pre-mRNAs (prp40-l) 
and 3' extended RNA species [rrp6A] (Supplemental 
Fig. 3). As expected, compared with the parental strain, 



the levels of intron-containing mRNAs — in particular, 
ribosomal protein mRNAs — were substantially reduced 
in the prp40-l strain (Supplemental Fig. 4A). This global 
reduction of spliced mRNAs in the prp40-l mutant was 
likely the result of degradation of pre-mRNAs that 
accumulated in this mutant (Bousquet-Antonelli et al. 
2000; Hilleren and Parker 2003). Despite this, the prp40-l 
data set contained a higher number of cDNAs that 
mapped to introns (Supplemental Fig. 4B), indicative of 
pre-mRNA accumulation due to defective splicing in 
this strain. 

Remarkably, in sharp contrast to intron-containing 
mRNAs, a large number of intronless mRNAs were up- 
regulated in the prp40-l mutant, indicative of a potential 
regulation by SMD (Supplemental Fig. 4A). However, we 
cannot rule out a possibility that increase in the levels of 
some of these intronless mRNAs in the prp40-l mutant is 
the result of pleio tropic effects. Compared with the parental 
strain, levels of 197 intronless transcripts increased at least 
threefold in the prp40-l mutant, and these were signifi- 
cantly enriched for genes involved in the unfolded protein 
response (P-value < 0.001). Twenty-four of 197 (-12%) 
contained predicted 5'ss and BPs in the right orientation, 
and 16 (—8%) were also enriched in the Sm-IP data 
(Supplemental Fig. 4C). Our analyses identified four 
transcripts [OYE3, SSA4, KAR2, and BDF2) that we pre- 
dicted were most likely to be SMD targets: They reproduc- 
ibly coprecipitated with the Sm proteins and contained 
splice signals in the correct orientation, and their steady- 
state levels increased at least threefold in the prp40-l 
mutant (RNA-seq and quantitative RT-PCR [qRT-PCR]) 
(Fig. 6A ; Supplemental Fig. 4C). 

Because SMD requires the activity of the nuclear RNA 
surveillance machinery, we looked for canonical splice 
junctions (GU/AG) in the sequencing data. The TopHat 
program (Kim and Salzberg 2011) identified splice junc- 
tions in the vast majority of known intronic mRNAs 
(79% [wild type] to 88% [rrp6A]) (Supplemental Fig. 5A). 
Splice junctions were also identified in 0.45% (wild type) 
to 1.1% (rrp6A) of mRNAs that are not known to contain 
introns. Based on these data, we predict that at least 1% 
of the intronless mRNAs might be subjected to SMD. For 
the vast majority of these transcripts, splicing would 
generate premature stop codons, and any spliced products 
would presumably be degraded. Indeed, in agreement 
with the unstable nature of SMD products, more splice 
junctions were identified in the rrp6A data set in addition 
to a higher number of reads mapped to protein-coding 
genes (Supplemental Fig. 5A). This indicates at least 
a partial stabilization of spliced products in the absence 
of Rrp6. Fewer splice junctions were recovered from the 
prp40-l data set, indicating that many of the identified 
junctions originated from a splicing event (Supplemental 
Fig. 5A). 

A total of 227 canonical splice junctions were found in 
49 transcripts, including SSA4, KAR2, and OYE3 (Supple- 
mental Table 2; Supplemental Fig. 5B), and 22 transcripts 
contained a BP consensus sequence within the predicted 
intron (Supplemental Table 3). Notably, splicing of OYE3 
at the identified splice sites would generate an in-frame 
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Figure 6. Other intronless mRNAs might be tar- 
geted by SMD. (A) KAR2, SSA4, and OYE3 transcript 
levels increase at least threefold in the prp40-l 
mutant. To validate the RNA-seq data, mRNA 
levels of the transcripts predicted to be targeted by 
SMD was measured by qRT-PCR. Error bars in- 
dicate standard deviations. (5) End point RT-PCR 
using oligonucleotides that span exon junctions 
confirms splicing of the KAR2, SSA4, and OYE3 
transcripts. (C) OYE3 spliced products accumulate 
in the rrp6A strain. Spliced and unspliced OYE3 
products were detected by Northern blot using 
a probe spanning the predicted intron. The BDF2 
exon 1 probe was used as a positive control. (D) 
Spliced products are stabilized in exosome (rrp6A) 
and (upflA) mutants. RNA from the RNA-seq ex- 
periments was subjected to qRT-PCR using oligo- 
nucleotides that span exon junctions. The data were 
normalized to the PPM2 reference gene. The Y-axis 
indicates the fold increase in the spliced RNA levels 
over the levels in the parental strain. Error bars 
indicate standard deviations. 
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deletion (Supplemental Fig. 6). However, as with BDF2, 
splicing could also occur at an alternative downstream 
3'ss, generating an out-of -frame product (Supplemental 
Fig. 6). Unfortunately, we did not find any BDF2 splice 
junctions in the RNA-seq data, but this could be due to the 
relatively low read coverage over BDF2. We did observe 
a modest increase in reads covering the predicted first exon 
in the rrp6A sequencing data (Supplemental Fig. 7), consis- 
tent with at least a partial stabilization of this product in the 
absence of Rrp6 (Fig. 2B). 

A subset of the genes contained putative introns in the 
5' untranslated region (UTR) (CPA2, GCY1, PIR3, SBH2, 
TBC3, YGR210C, and YR02) or 3' UTR (MSL5, RPL41A, 
QDR2, and YOR097C). Splicing of these mRNAs is 
therefore not expected to alter the reading frame (Sup- 
plemental Table 2) but could potentially affect mRNA 
stability and/or translation efficiency. 

Splicing of SSA4, OYE3, and KAR2 was confirmed by 
RT-PCR and sequencing of PCR products (Fig. 6B,D ; 
data not shown), and a strong accumulation of OYE3 
spliced products could be detected in the rrp6A strain by 
Northern blot (Fig. 6C). Interestingly, spliced mRNA 
products could be (partially) stabilized by deleting Rrp6 
or Upfl, indicating that, in these cases, both the nuclear 
exosome and cytoplasmic surveillance machineries con- 
tribute to their degradation (Fig. 6D). Because not all of 
the factors involved in cytoplasmic mRNA decay were 
analyzed here, it is possible that the actual percentage of 
intronless genes regulated by SMD is higher than we 
predicted (1%). 

We conclude that, in addition to BDF2, SMD also 
regulates expression of several other intronless genes. 



Discussion 

The Sm proteins play crucial roles in pre-mRNA splicing 
and exert their function by binding several spliceosomal 
snRNAs. We tandem affinity-purified RNAs associated 
with the Sm complex to identify novel Sm/splicesome 
targets. To our surprise, more than half of the top 500 
recovered mRNAs were not known to contain introns. 
Interactions between spliceosome components and non- 
intronic RNAs have been reported by other groups 
(Kotovic et al. 2003; Moore et al. 2006; Tardiff et al. 2006; 
Harigaya and Parker 2012). However, the biological 
relevance of these interactions remained unclear. Intrigu- 
ingly, recent studies have shown that the spliceosome is 
also used for purposes other than removing intron se- 
quences. In particular, it has been reported that fission 
yeast uses (L)Sm proteins and the spliceosome to mature 
the 3' end of the telomerase RNA, involving a single 
spliceosome-dependent cleavage step (Box et al. 2008; 
Tang et al. 2012). This encouraged us to investigate the 
association of the spliceosome with intronless transcripts 
in more detail. Our bioinformatics analyses show that 
many intronless mRNAs contain canonical splice signals 
in the correct orientation and therefore have the potential 
to recruit the splicing machinery. 

To investigate the potential function of the spliceo- 
some on intronless genes, we focused on BDF2, a gene 
encoding a bromodomain transcription factor, where we 
and others (Tardiff et al. 2006) observed recruitment of 
the spliceosome. Using yeast strains defective in nuclear 
RNA surveillance and pre-mRNA splicing, we discovered 
that splicing of the BDF2 mRNA led to an approximately 
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threefold reduction in BDF2 mRNA levels, and the resulting 
splicing products were rapidly degraded. This phenome- 
non, which we refer to as SMD, down-regulates mRNA 
expression as opposed to a well-documented positive 
outcome of conventional splicing and could potentially 
play an important role in down-regulating transcript 
levels. Indeed, the 5'ss and BP consensus sequences are 
conserved in genes homologous to BDF2 in related yeast 
species (Supplemental Fig. 8). Transcriptome sequencing 
of splicing and exosome mutants identified several splic- 
ing events on transcripts that were previously regarded 
not to have introns. We furthermore demonstrate that 
KAR2, SSA4, and OYE3 mRNA levels are also regulated 
via SMD and identify dozens of other potential targets. 
Finally, we show that Bdf 1, a Bdf2 paralog, is required for 
the recruitment of the spliceosome to BDF2, revealing a 
novel mechanism by which expression of paralogous 
genes can be regulated. 

A growing body of evidence suggests that defective or 
inefficient pre-mRNA splicing leads to the production 
of nonfunctional and potentially toxic RNAs such as 
unspliced mRNA precursors as well as various splicing 
intermediates. These are targeted by the cellular RNA 
surveillance systems for destruction. The current model 
implicates spliceosomal components in the initial iden- 
tification of defective splicing intermediates, which are 
discarded by the spliceosome and released to be degraded 
by the cellular RNA decay machineries (Houseley et al. 
2006; Egecioglu and Chanfreau 201 1 ; Parker 2012; Schmid 
and Jensen 2013). Unspliced RNAs are degraded by the 
exosome (from 3 ' to 5 ' ), Rat 1 /Xrn2, and Xrn 1 (from 5 ' to 3 ' ) 
in both nuclear and cytoplasmic compartments. Our 
results show that BDF2 splicing intermediates that 
accumulate due to delayed second step of splicing are 
also degraded. It was previously proposed that splicing 
intermediates are degraded in the cytoplasm (Hilleren and 



Parker 2003; Sayani et al. 2008; Mayas et al. 2010; Harigaya 
and Parker 2012). However, here we demonstrate that 
splicing products are also degraded by the nuclear exo- 
some, as they accumulate in cells lacking Rrp6. How are 
these RNAs recognized by the nuclear exosome? We 
speculate that recognition by the nuclear surveillance 
machinery is coupled to the nuclear retention times of 
these molecules. Growing evidence suggests that splicing 
and surveillance are tightly connected and that the fate of 
pre-mRNAs is dictated by competition between degrada- 
tion, splicing, and export machineries (Gudipati et al. 
2012; Sayani and Chanfreau 2012). Unspliced transcripts 
and splicing intermediates that are efficiently exported to 
the cytoplasm are presumably more likely to be degraded 
by the cytoplasmic RNA decay machinery, whereas 
transcripts with relatively slow export kinetics or long 
nuclear retention times are (also) targeted by the nuclear 
degradation pathway (Sayani and Chanfreau 2012). This 
model also implies that slow transition kinetics from the 
first to the second step of splicing could render RNA 
more susceptible to nuclear degradation. The efficiency 
of the 3'ss selection depends on several factors, includ- 
ing the distance of 3'ss from the BP, secondary structure 
of the intron, and interactions with the specific protein 
factors (Box et al. 2008; Meyer et al. 2011). Here we 
demonstrate that regulation of BDF2 by SMD can involve 
(at least) two different splicing events in which the same 
5'ss is used but two different BP sequences (see Fig. 7). We 
speculate that these two BPs are competing with each 
other for splicing factors and that selection of the BP is 
stochastic. If the spliceosome uses the upstream BP (+254), 
then it is very unlikely that the second splicing step will 
take place, as the BP is more than several hundred nucle- 
otides away from 3'ss elements (Zhang et al. 2007; Harigaya 
and Parker 2012). When the downstream BP is used, 
complete splicing can occur. Surprisingly, products gener- 
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Figure 7. Model for the regulation of BDF2 expression via SMD. In the presence of Bdfl, the spliceosome assembles on the BDF2 
mRNA. This leads to cleavage at the 5'ss and (1) subsequent release of exon 1 and (2) complete splicing. Splicing products are degraded by 
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ated by either pathway appear to be primarily degraded by 
the nuclear RNA surveillance machinery, as deletion of 
Rrp6 resulted in a substantial increase of partially and 
completely spliced products. Mutating either BP results 
in a twofold to threefold increase in BDF2 steady-state 
levels, suggesting that both pathways contribute signifi- 
cantly to SMD. The selection of BP and 3'ss signals is 
likely the rate-limiting step in splicing of BDF2, and any 
delay presumably leads to the activation of the nuclear 
surveillance machinery. This model is supported by a re- 
cent study from the Baumann laboratory (Kannan et al. 
2013) that focused on dissecting the mechanism of a one- 
step splicing reaction required for 3' end formation of the 
telomerase RNA in fission yeast. This study reports that 
slow kinetics of the second splicing reaction due to the 
suboptimal positioning of the 3'ss could lead to the 
structural rearrangements of the spliceosome complex 
and subsequent discard of the intermediates. We propose 
that, depending on the rate of the second step of splicing, 
the spliceosome could either direct RNA degradation [BDF2] 
or splice pre-mRNA (conventionally spliced mRNA- 
encoding genes). Importantly, each of these two scenarios 
leads to a completely different outcome: Production of 
unstable aberrant RNAs negatively impacts gene expres- 
sion (Fig. 7), whereas conventional splicing produces 
functional protein-coding mRNA. 

In S. cerevisiae, only a small proportion of genes 
contains introns and encodes transcripts that are spliced. 
As a possible explanation for this phenomenon, it was 
previously proposed that introns in many genes of the 
budding yeast were "lost" during the course of evolution 
(Fink 1987). However, recent studies as well as our data 
show recruitment of the spliceosome coinciding with the 
presence of conserved splice signals on genes that are not 
known to be spliced. It seems unlikely that these con- 
served sequences were retained through evolution with- 
out any functional implication. Indeed, we observed that 
a number of RNAs undergo splicing, creating fully spliced 
but nonfunctional RNAs that are subsequently degraded 
via SMD. While spliced BDF2 products are mainly targeted 
in the nucleus by SMD, some products of unproductive 
splicing can also be degraded by NMD in the cytoplasm. 

It is not clear exactly how many other genes are reg- 
ulated by SMD to the same degree as BDF2. Under the 
experimental conditions used, for some intronless tran- 
scripts in which we could detect splice junctions in the 
RNA-seq data, we could not detect a significant increase 
in pre-mRNA levels in splicing mutants, indicating that 
splicing of these mRNAs is inefficient. It is conceivable 
that the presence of (canonical) splice sites is sufficient to 
recruit the spliceosome and perhaps drive splicing to 
some degree, but for most transcripts, this is not suffi- 
cient to dramatically alter steady-state RNA levels. It is 
also possible that significant changes will only be detect- 
able under specific growth or stress conditions. Examples 
of this regulation were reported in higher eukaryotes. 
Thus, the well-documented unproductive outcome of 
alternative splicing can be used to regulate the expression 
of specific genes by suppressing the production of a pro- 
tein in the absence of the proper biological context 



(Lareau et al. 2007; Ni et al. 2007). For example, the 
expression of PTBP1 and PTBP2 factors is regulated via 
this mechanism in order to control their function during 
neuronal differentiation (Wollerton et al. 2004). In this 
scenario, alternative splicing results in production of 
a nonfunctional RNA degraded by NMD. Experiments 
in the laboratory are most frequently performed when 
cells are growing exponentially in high concentrations of 
glucose or galactose, conditions rarely found in nature. 
Under these conditions, cells devote the vast majority of 
their resources to the synthesis and splicing of pre- 
mRNAs encoding ribosomal proteins, masking the de- 
tection of less abundant or less efficient splicing events. 
Indeed, differential recruitment profiles of the spliceo- 
some and degradation machineries to target substrates 
were reported in response to different carbon sources 
(Bousquet-Antonelli et al. 2000), environmental stress 
(Moore et al. 2006; Pleiss et al. 2007), the cell cycle (Reis 
and Campbell 2007; Mullen and Marzluff 2008), and 
meiosis (Moore et al. 2006; McPheeters et al. 2009; 
Cremona et al. 2011). BDF2 expression has also been 
reported to be significantly affected under different stress 
conditions (Causton et al. 2001; Tkach et al. 2012). 

Consistent with recently published work, we also ob- 
served an increase in BDF2 mRNA levels upon deletion of 
another bromodomain-containing protein, BDF1, whose 
function in transcription initiation is redundant with 
Bdf2 (Matangkasombut et al. 2000; Durant and Pugh 
2007). We reasoned that the spliceosome could poten- 
tially regulate Bdf2 function in the cell by modulating its 
levels and therefore analyzed recruitment of Ul snRNP 
in bdfl A. Indeed, our data suggest that the approximately 
threefold increase in BDF2 mRNA levels (and approxi- 
mately twofold increase in protein level) in bdfl A is due 
to the reduced recruitment of the spliceosome to BDF2 
(Fig. 7). The function of Bdf 1 in regulating gene expression 
via splicing, in addition to its role in transcription initia- 
tion, may explain the multiple phenotypes associated with 
the lack of Bdfl. For example, BDF1 deletion causes slow 
growth, temperature sensitivity, salt hypersensitivity, 
flocculence in liquid culture, sensitivity to the DNA- 
damaging agent MMS, poor sporulation, and inability to 
grow on nonfermentable carbon sources (Lygerou et al. 
1994; Chua and Roeder 1995; Liu et al. 2007). 

It remains to be determined whether Bdfl is directly 
involved in the recruitment of the splicing machinery or 
whether this is a secondary effect of its function in 
transcription. Deletion of Bdfl was reported to affect 
snRNA levels, indicating a role in transcription of these 
genes (Lygerou et al. 1994); however, a more recent study 
has questioned these findings (Albulescu et al. 2012). We 
favor a model in which Bdfl plays a direct role in reg- 
ulating BDF2 expression via modulating spliceosome 
recruitment. The latter could be inferred from the fact 
that the presence of Bdfl on the BDF2 gene (Venters et al. 
201 1) has no noticeable impact on transcription (Fig. 4C). 
Moreover, loss of Bdfl was previously linked with ge- 
nome-wide defective splicing and impaired recruitment 
of the spliceosome to spliced genes, implying a general 
role in splicing (Albulescu et al. 2012). As regulation at 
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the level of chromatin has recently emerged as an additional 
regulatory mechanism affecting splicing (de Almeida and 
Carmo-Fonseca 2012), it would be interesting to examine 
whether Bdfl function in splicing is mediated via recog- 
nition of the acetylated his tones. 

In conclusion, we propose that the cooperative contri- 
bution of multiple mechanisms during RNA proces- 
sing — such as recruitment of the spliceosome, the check- 
point allowing the second step of splicing, and finally, 
recruitment of the exosome to the defective RNAs — is 
important to ensure dynamic regulation of gene expres- 
sion in eukaryotes. 

Materials and methods 

Yeast strains 

The S. cerevisiae strains used in this study are listed in 
Supplemental Table 4. A list of the oligonucleotides used is 
given in Supplemental Table 5. 

Northern blotting 

RNA extractions and Northern blot experiments were performed 
as previously described (Vasiljeva and Buratowski 2006). Gene- 
specific PCR-generated fragments were used as probes using the 
oligonucleotides listed in Supplemental Table 5. 

Purification of the Sm complex 

TAP-tagged Sm complex was purified from 16 L of yeast culture 
essentially as described in Vasiljeva and Buratowski (2006). RNA 
was extracted from calmodulin -bound material and subjected to 
RNA-seq. 

Transcriptional run-on 

Transcriptional run-on was performed as described in Birse et al. 
(1997). Single- stranded probes for transcriptional run-on analysis 
were generated by in vitro transcription using T7 or SP6 poly- 
merase from plasmids harboring regions corresponding to differ- 
ent regions of BDF2. BDF2- and ADH1 -derived fragments were 
generated using the primers listed in Supplemental Table 5 and 
were cloned into the pCR-Blunt II-TOPO vector (Invitrogen). 

ChIP experiments 

ChIP procedures and quantification were performed as described 
(Keogh and Buratowski 2004; Kim et al. 2004). For temperature 
shift experiments, cells were incubated at 23°C until OD(A 60 o) = 
0.5, and cells were further incubated for 120 min at 37°C. The 
primers are listed in Supplemental Table 5. 

High-throughput RNA-seq 

Procedures for cDNA library preparation and RNA-seq are de- 
scribed in the Supplemental Material. Raw (fastq) and processed 
(bedgraph and GTF files) RNA-seq data can be downloaded 
from the NCBI Gene Expression Omnibus repository (http:// 
www.ncbi.nlm.nih.gov/geo, accession number GSE49966). 
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